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We introduce a topological flux function to quantify the topology of magnetic braids: non-zero, line-tied 
magnetic fields whose field lines all connect between two boundaries. This scalar function is an ideal invariant 
defined on a cross-section of the magnetic field, and measures the average poloidal magnetic flux around any 
given field line. Moreover, its integral over the cross-section yields the relative magnetic helicity. Using the 
fact that the flux function is also an action in the Hamiltonian formulation of the field line equations, we 
prove that it uniquely characterizes the field line mapping and hence the magnetic topology. 



I. INTRODUCTION 

In many plasmas, ranging from astrophysics to mag- 
netic confinement fusion, the topology — i.e., the linking 
and connectivity of the magnetic field finest — is an ap- 
proximate invariant of the dynamics. This is because 
these plasmas typically have such low dissipation that, to 
first approximation, their evolution is ideal. That is, on 
large scales where the magnetohydrodynamic approxima- 
tion holds, they satisfy an ideal Ohm's law-2, preserving 
the magnetic topology. Therefore a practical question is, 
given two magnetic fields satisfying the same conditions 
on the boundary of some volume V, can one be reached 
from the other by some ideal evolution in V? 

We restrict our attention to line-tied magnetic flux 
tubes, where all field lines stretch between two bound- 
aries and the magnetic field in the volume is non- 
vanishing. This models, for example, a coronal loop 
in the Sun's atmosphere, where the footpoints re- 
main essentially fixed on the rapid timescale of coronal 
relaxation*. To simplify the presentation in this paper, 
we consider a magnetic field B defined on a cylinder 
V = {(r,(f>,z)\0 < r < R,0 < z < 1}, satisfying B z > 
everywhere in V, and impose the boundary conditions 
that B|ay = e z and v\gy = 0, where v is the veloc- 
ity. Extensions of the results to more general boundary 
conditions are discussed in Section [Vj For convenience, 
we call magnetic fields satisfying the above conditions 
"magnetic braids" (Fig. 1). Two magnetic braids are 
topologically equivalent if they can be linked by an ideal 
evolution where v = on dV throughout. 

In principle, one can determine whether two magnetic 
braids are topologically equivalent by comparing their 
field line mappings from z = to z = 1. Throughout 
this paper, let f(xo; z) G V denote the point at height z 
on the field line traced from xq = (r , (f>o, 0) on the z = 
boundary. Under this parameterization of field lines by 
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z, we have 

df(x ;z) = B(f(x ;z)) 
dz B z (f(x ;z))' 

For shorthand, we shall denote the mapping from z = 
to z = 1 as F(xq) = /(xq; 1). Under our boundary con- 
ditions, two magnetic braids B, B are equivalent if and 
only if F = F. Note that if we were to relax the condi- 
tion that B = e z on the side boundary r = R, then F 
would determine the topology only up to an overall rigid 
rotation through 2mr, n£Z. Symplectic mappings have 
long been used themselves as models of periodic magnetic 
fields in fusion devices^, and field line mappings have also 
been used extensively for characterizing line-tied coronal 
magnetic fields^. But the mapping is usually very sen- 
sitive to small fluctuations in the underlying magnetic 
field. This makes it very difficult, if not impossible, to 
determine whether two field line mappings can be related 
by an ideal evolution in anything other than highly ide- 
alized situations. 

A much more robust topological quantity is the total 
magnetic helicity, which has a broad range of applications 
in both laboratory and astrophysical plasmas^. It has 
proved so robust that it has been hypothesized to be 
the only quantity determining the final state of turbulent 
relaxation in reverscd-ficld pinches and similar devices^. 
But the helicity is an extreme reduction of the topological 
information in the three-dimensional magnetic field to a 
single number, and it does not uniquely characterize the 
topology. There are a large class of magnetic fields that 
have the same helicity but different field line mappings. 

In this paper, we describe a quantitative measure 
A(xo) which we call the "topological flux function" . This 
is a scalar function defined on a cross section of the mag- 
netic field, or cquivalently on each field line. It is more 
robust than the field line mapping, while containing more 
detailed information than the helicity. A similar func- 
tion was introduced for magnetic fields in a half-space 
by Berger— , who showed that it is effectively a helic- 
ity per field line. The function A has also appeared in 
the literature in a different guise: as an action integral 
yielding the magnetic field line equations in a variational 
formalism^^. These interpretations arc discussed in 
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more detail in Section Mil Recently, we have described 
how A may also be viewed as a generalization of the 
scalar flux function used to define two-dimensional mag- 
netic fields, and how it may be used to define and mea- 
sure magnetic reconnectionii. Here we go further. Our 
main result is that, with a particular choice in its defini- 
tion, A uniquely characterizes the topology of a magnetic 
braid. In other words, A = A for two magnetic braids 
if and only if they are topologically equivalent. More- 
over, not only can A determine whether two braids are 
topologically equivalent, it can also quantify how much 
dissipation or reconnection is needed to connect the two 
states. This will be invaluable in future dynamical stud- 
ies of such systems. 

This paper is organized as follows. In Section[Hl we de- 
fine A and give its basic physical interpretation in terms 
of magnetic flux. Its interpretations as a field line helic- 
ity and as a Hamiltonian action are described in Section 
Mil In Section HVl the Hamiltonian interpretation is used 
to prove our main result that A uniquely characterizes 
the magnetic topology. We outline in Section [V] how the 
boundary conditions may be relaxed. 



II. DEFINITION OF THE TOPOLOGICAL FLUX 
FUNCTION 

The topological flux function A is defined simply by 
integrating the vector potential A (where B = V x A) 
along magnetic field lines. It may be written as 



A(x ) = 



■=i 



f' A(/(xo;z)).B(/(xo;z)) 
B z (f(x ;z)) 



(2) 

Broadly, A is conceived to measure poloidal (horizon- 
tal) magnetic fluxes in the domain V. But it will be 
beneficial to restrict the gauge of the vector potential A 
in the definition @. This is because the value A{xq) for 
a general field line is not gauge invariant. Indeed, under 
a gauge transformation A — > A + Vx, 



A^A + F* X -X, 



(3) 



where we use the pull-back notation F*x(xq) = 
x(F(xo))- Provided that x is chosen to be periodic 
in z, then it follows that for periodic field lines, where 
F(xq) = xq, the value of A(xq) is gauge invariant. In an 
ideal evolution, two periodic field lines define a comoving 
surface with flux A(xi) — A(xq), as in Fig. []Ja), and the 
set of periodic field lines can be used to define a poloidal 
flux partition of the magnetic fieldii. But for a general, 
non-periodic field line, A is not gauge invariant unless 
we impose further gauge conditions on A in its definition 
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We show here that A is a physically meaningful quan- 
tity for any field line if we impose 
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FIG. 1. (Color online) Magnetic braids in the cylinder V, 
showing (a) the comoving surface defined by two periodic field 
lines, and (b) the poloidal surface at viewing angle <j> for a 
single field line (bounded by the arrowed path). 



where A = (r/2)e r p is the vector potential of a refer- 
ence field B rcf = e z that matches B on dV. A similar 
restriction is used to ensure gauge invariance in the well- 
known relative magnetic helicity-i^. To demonstrate how 
A becomes meaningful, define a "poloidal surface" with 
flux $(0) bounded by (i) the field line, (ii) a vertical line 
on the side boundary r = R at azimuth <f>, (iii) a straight 
line Lq on z = joining the startpoints of the first two 
lines, and (iv) a straight line L\ on z = I joining their 
endpoints (Fig. [lb). In view of the boundary condition 
v\gy = 0, this is a comoving surface and $(</>) is an ideal 
invariant for any </>. By Stokes' Theorem, 



$(0) = A(x ) 



A - dl, 



L +L x 



(5) 



with no contribution from the side boundary since A z = 
there by our gauge choice. For a periodic field line, 
the integrals in ([3]) along Lq and L\ will be equal and 
opposite, so the flux $((/>) is independent of the "viewing 
angle" (f> and is given by A(xq). For a non-periodic field 
line, the integrals need not cancel, and will depend in 
general on <f>. However, we can show that in this case 
A(xq) gives the average flux over all viewing angles. To 
see this, consider the shaded quadrilateral lying on the 
lower boundary in Fig. [Tf b) . The integral over Lq and 
L\ returns the magnetic flux through this quadrilateral, 
since A T = on z = 0, 1. Since B z = 1, this is simply 
the area of the quadrilateral, 

/ A- d\= -{i^sin(F -cf>) + r o sin((/>o-0)}. (6) 

This expression vanishes on averaging over <f> so, for this 
gauge restriction, 



1 

2^ 



$(</>) d<j> = A(x ) 



(7) 



n x A 



dV 



n x A 



rcf 



\dV, 



(4) for any field line. 



A fundamental property of A in our restricted gauge 
((3]) is ideal invariance. One can see this by calculating 
dA/dt as a line integral over a moving domain^. In an 
ideal evolution, dA/dt = v x B + Vip, so 



dA 
~dt 



z=0 
z=l 



dA 
~dt 



vw>- 



-v x V x A + V(v- A) | • dl (8) 
v • A) • dl. (9) 



The second term vanishes since ~v\dv = 0. Under our 
gauge restriction (j4j , ip is constant on dV and the integral 
vanishes. 



III. PHYSICAL INTERPRETATIONS 

Before proving our main result, we describe two illu- 
minating interpretations of the topological flux function. 



A. Field Line Helicity 

There is a simple relation between A and the magnetic 
helicity^. Since V is magnetically open, we use the rela- 
tive helicity H r (Ref. Il2l ). Under our gauge conditions, 
however, this reduces to the same expression as the total 
helicity, 

H r = [ (A + A rcf )-(B-B rcf )d 3 x = I A-Bd 3 x. (10) 
Jv Jv 

Now suppose we change variables from (r, <f>, z) to 
(r ,4> ,z), where (r,(j>,z) = f(x ;z) and a- = (ro,0o,O) 
is the footpoint on z — of the field line through (r, <fi, z). 
The Jacobian of this coordinate transformation is 



det(J) 



r a B z (x ) 
rB z (r 7 <f>,z) 



(11) 



as may be verified by considering a thin flux tube around 
the field line and using V • B = 0. Thus we can re- write 



H r = I A(r,(p, z) ■ B(r, (f>, z) rdrd<fidz 
'v 



(12) 



A{f(xo;z)) ■ B(f(x ;z)) - f z ^ X °^ N , d 2 x c 



B z (f(x ;z)) 



A{x )B z (x )d 2 x . 



(13) 
(14) 



z=0 



So A is a density for H r in the cross-sectional plane, 
weighted by magnetic flux. (With our boundary condi- 
tions, we simply have B z (xq) = 1.) Hence Berger— calls 
A a "field line helicity" . Clearly it is possible for A / 
even when H r = 0, providing that the integral of A over 
all field lines vanishes. 

An example of this is the magnctohydrodynamic sim- 
ulation of a relaxing solar coronal loop described by 




FIG. 2. (Color online) Example calculation of the topological 
flux function A in a numerical magnetohydrodynamic simu- 
lation of magnetic relaxations^, using the gauge conditions 
and A r = 0. Color contours on the lower boundary show 
A. Panels (a) and (b) correspond to before and after the 
relaxation respectively. 



Wilmot-Smith et aZ.— ^ and shown in Fig. [2] Although 
the evolution is resistive and there is widespread recon- 
nection as the field relaxes, the relaxation is sufficiently 
fast to preserve the initial helicity H r = 0, as predicted 
by Taylor theory^. But contrary to Taylor theory, which 
would predict a uniform relaxed field B = e z with .4 = 0, 
the final state maintains equal regions of positive and 
negative A, manifesting itself in non-trivial topology of 
the field lines in the end-state of the relaxation, de- 
spite the conservation of helicity^ (see also Ref. liH ). 
This illustrates how A contains more detailed informa- 
tion about the topology than H r . 



B. Hamiltonian Action 

To motivate why A might be a sufficient condition 
for two magnetic braids to be topologically equivalent, 
we point out another physical interpretation of A that 
demonstrates a deep connection with the magnetic field 
structure. Namely, A is the action in a variational for- 
mulation that leads to the equations of the magnetic field 
lines^. In other words, for given vector potential A and 
field line mapping, the magnetic field lines x(Z) are given 
by the Euler-Lagrange equations that extremize the in- 
tegral 



^ A(x).— dl. 



(15) 
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It is well known that the magnetic field lines arc tra- 
jectories of a Hamiltonian system — a fact well exploited 
in the modeling of toroidal fusion devices^. So the ac- 
tion (fT5|) is that of a Hamiltonian system, although it 
is not necessarily written in canonical coordinates^. To 
demonstrate that it really is Hamiltonian, we can re- 
write (|15j) in canonical coordinates by a gauge transfor- 
mation A — > A + V%. This adds dx/dl to the integrand 
(the Lagrangian), which leaves the field line equations 
(the Eulcr-Lagrangc equations) unchanged. We can then 
write our system in canonical coordinates by choosing an 
appropriate gauge. If we set A r = everywhere in space, 
then A • dl = rd(j> + A z dz. Making the identifications 
po r^(r,^z), <7 ^ </>, t «-> z, H o —A z (r, <j>, z), 
our action becomes 



A = 



This is a 1 degree-of-frcedom Hamiltonian system in 
canonical form. The generalized coordinate is <p, the gen- 
eralized momentum is rA^, and the Hamiltonian is — A z . 
Time corresponds to our z coordinate, so our Hamilto- 
nian is in general time dependent. We remark that the 
canonical gauge choice A r = is equivalent to writing B 
in the form B = V(rA^) x V</> + V(A Z ) x Vz, which is 
widely used in toroidal plasmas^. This gauge is also con- 
sistent with the gauge restriction Q imposed to ensure 
ideal invariancc of A. 



IV. UNIQUE TOPOLOGICAL CHARACTERIZATION 

Having identified A as the action in a Hamiltonian 
system, we now use general results about Hamiltonian 
systems to prove our main result that A is not only a 
necessary condition but also a sufficient condition for 
two magnetic braids to be topologically equivalent. The 
argument is based on work of Haro^ on primitive func- 
tions of exact symplectomorphisms, and is best expressed 
in terms of differential forms. We assume that A is in 
canonical gauge A r — and satisfies 

As a special case of a more general result about Hamil- 
tonian systems, it follows that 



dA = F*a - a, 



(17) 



where a = (r 2 /2)d(f> is the canonical (Liouville) 1-form 
(Ref. Llll p. 148). For more details and proof of (fTT]) . see, 
for example, Ref. HO (Proposition 9.18) or Haro^, where 
A is called a "primitive function" for F. Writing out the 
two components explicitly, (|17[) says that 



dA 

dr 



(F r ) 2 \ OF, 



dA 



dr d(f> 



(F r ) 2 \ 3F t 



d(f> 



(18) 

To prove that A uniquely determines the topology, sup- 
pose that we have two magnetic braids B, B, with topo- 
logical flux functions A, A and field line mappings F, 



F respectively. We already know that F = F implies 
A = A, since A is an ideal invariant. But we can also sec 
this from (|17p . which gives 



dA = F*a 



F*a 



a 



dA, 



(19) 



so that A and A differ by at most an overall constant, 
which vanishes since both braids satisfy A(R, <f>) = 0. 

To prove the converse, assume that A — A, and define 
the mapping G = F o Then, using (fTT). 

G*a — a= (F -1 )* o F*a — a (20) 

= (F- 1 )*{a + dA)-a (21) 

= (F^Ya-a+iF^ydA (22) 
= (F~ 1 )*a - (F- 1 )* o F*a + (F'^dA (23) 

= (F- 1 )*(a-F*a) + (F- 1 )*dA (24) 

= -(F~ 1 )*dA + (F _1 )*cM (25) 

= 0. (26) 

Now we determine the possible mappings G satisfying 
G*a — a, or equivalently (compare Eq. ([T8"]0 



(G r ) 2 dGj 

2 dr 



dr 



(G r ) 2 OCp r\ 



#o = 0. (27) 



The tq and <f>o components give, respectively, G$ = g(4>o) 
and G r = r$(dg / 'dip^)" 1 / 2 . But the possibilities arc re- 
stricted by our boundary conditions: firstly, G r (R, <j>o) = 
R implies that dg/depo = 1, so that G is a rigid rotation 
G(ro,0o) = (r ,(f>o + const). But then G<f,(R,4>o) = 4>o 
implies that G = id, and so F = F. This completes the 
proof. 



V. DISCUSSION 

For clarity of presentation, we have assumed certain 
boundary conditions on the magnetic field, namely that 
B = 1 on all boundaries of our cylinder V, and that 
v = on the bottom and top boundaries. We indicate 
here how the results generalize when these conditions are 
relaxed. 

If one allows ^ on the side boundary r = R, 
then G is determined only up to an overall rigid rotation, 
and we would have the result that A and A differ by a 
constant if and only if F and F differ by an overall rigid 
rotation. Such a rotation can be detected from knowledge 
of F and F on the side boundary alone. In physical 
applications, the mapping on the side boundary may well 
be fixed in time. 

It is also possible to consider more general B z distri- 
butions on the boundaries z = 0, 1, providing A satisfies 
(|4]) for an appropriate A ref . In that case, the canonical 
1-form a is a = rA ( j ) (r,(f>,z)d<f), and the possible map- 
pings G that preserve a (and are undetected by ^4) may 
differ, although they take the general form of a cotangent 
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lift G(<f) ,p) = (T((f>o),p(dT/d(j)o)) 1 wheiep = rA^ is the 
generalized momentum (see Ref. [l9|, Proposition 6.3.2). 
If F = F on the side boundary r = R, it follows that 
T(</>o) = 4>o and hence again G = id. 

Finally, one could consider a toroidal domain where B 
and v are periodic in z and remove the condition that 
v = on the boundaries z = 0,1. In that case, the 
freedom to shuffle around field line endpoints on z = 0, 1 
means that topological equivalence is a weaker notion, al- 
though it is certainly not true that any two field line map- 
pings are equivalent, as would be the case if one had the 
freedom to apply independent motions on both bound- 
aries. In the periodic case, the cross-section z = 0, 1 
is no longer distinguished by the boundary conditions. 
Changing cross-section has the same effect on F as an 
ideal evolution, so that the new mapping may be written 
F = S o F o S^ 1 for some field line mapping S. If A is 
the flux function for F, then one can show that the flux 
function A for F is given by 

A=(S- 1 )*(A + F*x~x), (28) 

where y is related to the mapping S by S*a — a = d\ (see 
Ref. Il81 ) . The practical problem of determining whether 
two given topological flux functions are related in this 
way remains for further investigation. 
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